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SUMMARY 

An improved method of analyzing multicomponent rotating assemblies for 
the determination of operating stresses and component compatibility has been 
developed. In this method, a single finite element model is developed which 
contains all of the separate components in the rotating assembly. This is 
made possible by using gap elements to simulate the contact surfaces between 
components. The MARC finite element computer program is then used to perform 
the analysis. 

This improved method is less time consuming and more reliable than the 
conventional method of analyzing such systems. Previously, a method called the 
flexibility method was used at the NASA Lewis Research Center for this task. 

It requires the individual analysis of components, derivation of compatibility 
equations, determination of certain coefficients for these equations, and an 
iterative solution procedure. 

This report presents results for two different stress-compatibility analy- 
ses of a six-component axial flow compressor rotor. The results for the flex- 
ibility analysis method are compared with those for the improved analysis 
method. It was found that the stresses predicted by each method compare quite 
well with each other. The predictions of the component compatibility, as well 
as the magnitude of the forces at the contact surfaces, also compare well for 
these two analysis procedures. It is therefore recommended that the improved 
analysis method be used to determine the stress-compatibility characteristics 
of multicomponent rotating systems. 


INTRODUCTION 

Much of the rotating hardware in turbomachinery - such as the rotors of 
compressors, turbines, pumps, and blowers - is designed as an assembly of sev- 
eral components which are mated together by means of an external-internal pilot 
interface. A typical pilot design is shown in figure 1. A tight pilot inter- 
face is essential in rotating machinery since it maintains the concentricity 
of separate components that are rotating about the same axis. In some cases 
the engineer may take advantage of the supporting characteristics of the pilot 
interface - the outer (enclosing) part provides radial support to the inner 
part. However, as an assembly of components is subjected to certain loading 
conditions, such as centrifugal loading, one of the outer parts may expand (in 
the radial direction) more than its adjacent inner part. This can loosen the 
parts and upset the integrity of the assembly. One of the difficulties in 
designing a pilot feature is that the pilot surfaces must be orientated such 
that the inner part expands at a greater rate than the outer part. In some 





cases, a sufficient assembly interference fit can be chosen that ensures that 
the parts will not separate when subjected to centrifugal loads. 

The rotating assembly of parts is usually held together by axial pre- 
loaded bolts or rods. Multiple tie bolts are sometimes used in a circular pat- 
tern. A single tie bolt, which lies along the rotational axis of the machine, 
can also be used. 

Rotating hardware that is designed in such a manner is subjected to sev- 
eral loading conditions while in operation. These conditions include centrifu- 
gal forces due to rotation, axial forces due to the preloaded bolts or rods, 
and pilot surface forces due to initial pilot interference fits. Other load- 
ing may include forces caused by transient or steady-state thermal distribu- 
tion, pressure forces, or aerodynamic forces. 

To properly determine the stresses and deflections of this type of assem- 
bly, the analyst cannot simply analyze each component independently and insure 
that the pilot interfaces are designed in the proper external-internal orienta- 
tion. The forces at each contacting surface are dependent on the forces at any 
or all of the other contacting surfaces within the assembly. For example, the 
increasing pressure at one pilot interface (due to centrifugal loading) may 
cause the pressure at an adjacent pilot interface to decrease, possibly forming 
a gap between two parts. Another consideration is that the pilot pressures, as 
well as the axial clamping force due to the tie bolt preload, contribute to the 
stress distribution on the assembly. 

The design and analysis are complicated further because several considera- 
tions must be taken into account when selecting the tie bolt preload. For 
instance, as an assembly of rotating parts increase in speed, the Poisson's 
ratio effect tends to shrink the axial dimensions of the assembly, thus reduc- 
ing the load on the bolt holding the components together. If the preload in 
the bolt diminishes to zero, the parts will loosen, causing an unbalance and 
possibly catastrophic results. Therefore it is very important to be able to 
accurately determine what preload remains in the bolt at operating conditions. 

The conventional method of analysis that has been used at NASA Lewis for 
this type of assembly is based on the use of the principle of complementary 
energy. This conventional method requires that each component of the rotating 
assembly be analyzed separately (usually a finite element analysis) under a 
variety of loading conditions in order to determine the interdependence of the 
contacting surfaces, such as the displacement of one surface due to a load at 
another surface. These interdependencies are described as influence coeffi- 
cients, and are determined for each component in the assembly. An equilibrium 
equation is then derived for every contacting point in the assembly, by equat- 
ing the displacements of the two surfaces in contact. These equations contain 
the previously determined influence coefficients, displacement at the contact 
surfaces due to operating loads, interference fits of the pilots, and pressure 
forces at the interface of the contact points. The interface pressure forces 
are the unknowns in each equation and must be determined by solving a system 
of N equations and N unknowns where N is the number of contacting points 
in the assembly. Once the system of equations is solved, the analyst can 
determine if contacting surfaces remain closed. If gaps develop, or if the 
tie bolt cannot hold the assembly of components together at operating condi- 
tions, then design modifications must be made and the analysis repeated. Aside 
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from being very time consuming, this type of analysis involves the manipulation 
of large quantities of numbers and is very prone to mistakes. Errors are not 
obvious and may go unnoticed throughout the analysis, giving false results. 

Other methods for analyzing this type of assembly exist (ref. 1), and they 
usually utilize the principle of minimum total potential energy. For all of 
the methods, a set of equilibrium equations must be solved to determine the 
unknown forces of a statically indeterminate structure. 

An improved method of analyzing multicomponent rotating systems has been 
developed which utilizes the MARC finite element computer program. By using 
MARC gap elements to simulate the points of contact, the entire assembly of 
components can be represented in a single finite element model. Therefore it 
is not required to develop and solve a system of equilibrium equations, nor to 
analyze each component separately for the determination of influence coeffi- 
cients. Also, any changes that must be made in the design are simply done by 
modifying the model that represents the complete assembly. This method signif- 
icantly reduces the time and effort required for such an analysis. 

The objective of this report is to describe the improved method of analy- 
sis. A stress-compatibility analysis of a seven-component axial flow compres- 
sor rotor has been performed, using the two methods described above, and the 
results are compared. 


SYMBOLS 


du incremental displacement vector 

df incremental force vector 


K 


stiffness matrix 



N 


displacement at point i in the direction of Pi due to a 
load Pj (influence coefficients)^ 

number of contacting points in the assembly 


Pf 3 force due to the tie bolt 

Pi contact force at interface i (integrated around the circum- 

ference) 


Pl.P3.P5.P7.P9 radial pressure forces corresponding to the pilot interfaces 

P2.P4.P6.P8 forces that tend to rotate the radial contact surfaces 
between components 


51 shrink (in axial direction) of all rotor components, due to 
radial pilot forces and compression from the tie bolt 

52 shrink (in axial direction over the length of the tie bolt) 
of all rotor components, due to centrifugal and thermal 
loading 
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S3 shrink (in axial direction) of the tie bolt, due to centrifu- 

gal and thermal loading 

AB stretch of the tie bolt caused by the axial force, Pg, in the 

tie bolt 

AL total initial axial stretch in the tie bolt that would be 

required to cause the axial forces, Pg, to be just equal to 
zero 

Ayi initial pilot interference fit in the direction of P^ 

Si displacement at point i in the direction of P^ due to cen- 

trifugal and thermal loading onlyl 

Superscripts: 

D1 , D2,D3,D4 disks 1 , 2, 3, 4 

FS front shaft 

RS rear shaft 

TB tie bolt 

Subscripts: 

i,j designations for i and j are 1 to 9 and B 


DESCRIPTION OF COMPRESSOR 

The compressor rotor that was analyzed is from the three-stage axial flow 
compressor shown in figure 2. This is an aerodynamic research compressor 
which is driven by a variable-speed electric motor through a speed increasing 
gearbox. The rotor design includes a single tie bolt and is composed of six 
separate components: a front shaft, four disks, and a rear shaft (fig. 3). 

Compressor blades are integrally machined into the first three disks, and the 
fourth disk is used as an axial spacer. The blades were not included in this 
analysis, but were left for a separate detailed analysis, although the centrif- 
ugal forces caused by them were considered. 

The four disks are mated together by means of a external-internal pilot 
interface, and the pilot diameters are dimensioned such that there is a 
0.0254-mm (0.001-in.) radial interference fit between the mating parts. The 
disks are axially bolted to the rear shaft by a single tie bolt which lies 
along the axis of rotation. During assembly, the tie bolt must be torqued 
such that its preload will prevent gaps from forming at the axial contact sur- 
faces between mating parts. After the tie bolt is installed, the front 


displacements due to forces P 2 , P 4 , P 6 » and Pg are taken to be the 
differences in the axial displacements at the outer and inner points where the 
forces are applied. 
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shaft is mated to the first disk by means of a pilot, and it is axially bolted 
with four bolts located on a 63.5 mm <2.5 in.) diameter. 

The compressor components are constructed of the following materials: 
titanium alloy T i - 6 A 1 -4V for the four disks, AISI 4340 alloy steel for the 
front and rear shafts, and 18 percent nickel alloy maraging steel (200 ksi) 
for the tie bolt. The material properties that were used in the analysis are 
shown in table I. 

The compressor rotor was analyzed at an operating speed of 47 500 rpm. 

The airflow through the compressor varied in temperature from 15 °C (60 °F) at 
the inlet to 182 °C (360 °F) after the third stage of compressor blades. The 
operating temperature for the bearings was 82 °C (180 °F). 


FLEXIBILITY ANALYSIS METHOD 

This particular analysis technique is sometimes referred to as the flexi- 
bility method and is based on the principle of complementary energy. Refer- 
ences 2 to 4 describe the theoretical background pertaining to the principle 
of complimentary energy, and a general description for the applications of this 
method is given in references 5 and 6 . The analysis involves the formulation 
and solution of a set of equilibrium equations hereafter referred to as compat- 
ibility equations. The general equations are derived using the principle of 
superposition in which the deflections of any two components, due to all 
applied known and unknown loading, are equated at their interface. One equa- 
tion is derived for each contacting surface in the assembly of components that 
make up the compressor rotor. The satisfaction of these equations indicates 
that the surfaces are compatible and that no gaps develop between parts. The 
equations are written in matrix form, and a simultaneous solution must be per- 
formed to solve for the unknown forces. The solution may indicate impossible 
loads, such as a tensile force across a contacting surface. In this case, 
the tensile forces must be set equal to zero, and another solution must be 
obtained. This iterative process is complicated by the fact that there is an 
interdependence of the effects of the unknown forces. 

The FEATS finite element computer program (ref. 7) was used to analyze 
the seven components of the compressor rotor separately in order to determine 
the values of certain coefficients contained within the set of compatibility 
equations. Two-dimensional four-noded quadrilateral axi symmetric elements 
were used in the modeling. The finite element models, along with the unknown 
interacting forces, are shown in figure 4. Since there are 10 unknown forces 
on the assembly, the compatibility equations will be described with a 10 by 10 
matrix. The forces P] , P3, P 5 , P7, and Pg are radial pressure forces corre- 
sponding to the pilot interfaces. Forces that tend to rotate the radial con- 
tact surfaces between components are described as P2, P4, P6, and Pg. Notice 
that these forces only create a moment on the surfaces. The force that is 
imposed on the components due to the tie bolt is Pg. 

By equating the displacements at each of the contacting surfaces (1 to 9), 
the following compatibility equations were derived: 

Surface 1 : 


K 11 P 1 ♦ 5 1 S ■ K n p i " K ?2 P 2 ♦ k u p 3 + k ?b p b + 


,D 1 




( 1 ) 
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Surface 2: 


P + P 
^21 P 1 K 22 2 


+ K 23 P 3 + K 2B P B + & 2 


+ K^P + K^P + K^P + + £^2 

^22 P 2 R 23 H 3 K 24 P 4 ^25 P 5 K 2B P B 6 2 


( 2 ) 


Surface 3: 

K^^P + P 
*31 1 32 P 2 


* K ?3 P 3 * K 5b P B + S “' 


l^D2 d ^D2 d ^D2 d ^D2 d „D2 d -D2 Aw 

K 32 P 2 K 33 P 3 + K 34 P 4 + K 35 P 5 + K 3B P B + S 3 + Ay : 


(3) 


Surface 4: 

K^P + K^P 
K 42 P 2 43 3 


^D2 d ^D2 d ^D2 d -D2 
+ K 44 P 4 + K 45 P 5 + K 4B P B + S 4 


K^^P + K^P + K^^P + K^^P + K^P + £^ 
^44 P 4 + K 45 P 5 K 46 P 6 + K 47 P 7 + K 4B P B + 6 4 


(4) 


Surface 5: 

Kf 2 ? 2 k“p 3 * 


K 54 P 4 + K 56 P 6 + K ?B P B + 


l/D3q , u-D3p ^D3 d ^D3 d ^D3 d c D3 a 

K 54 P 4 K 55 P 5 K 56 P 6 + K 57 P 7 + K 5B P B + 5 5 + Ay 5 (5> 


Surface 6: 

K^^P + K^P 
K 64 P 4 + ^65 P 5 


+ K 66 P 6 + K 67 P 7 + K 6B P B * 5 f 


i(D4p i/D4p t/D4p + K^^P + K^^P + 

R 66 P 6 o7 P 7 K 68 P 8 K 69 P 9 K 6B P B 6 6 


(6) 


Surface 7: 

K^P + K^P 
K 74 P 4 + *75 P 5 


+ K ?6 P 6 + K ?I P 7 * K ?B P B + 5 “ 


^D4 d ^D4 d .,D4 d „D4 d ^D4 d r D4 . 

K 76 P 6 + K 77 P 7 + K 78 P 8 + K 79 P 9 + K 7B P B + 6 7 + Ay 7 


(7) 


Surface 8: 

K^P + K^P 
K 86 P 6 87 P 7 


, K^^P + K^P + 

+ l^ggPg + (S gg P g + kgg‘ g + 6g 


l^RSp i/RSp RS p pRS 

^88 P 8 ^89 P 9 aB P B + 6 8 


( 8 ) 
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Surface 9: 

K^p + + K^P + K^P + K^P + 

^96 H 6 ^97 P 7 R 98 M 8 ^99 K 9 Sb P B 6 9 


.RS D .RS d .RS d c RS a 

K 98 P 8 + K 99 P 9 + K 9B P B + S 9 + Ay 9 


( 9 ) 


The compatibility equation for the contact surface located at the head of 
the tie bolt was derived as follows: 


AB = AL - SI - S2 + S3 


( 10 ) 


where 


k b! p i + 

C P 2 + 

K B3 P 3 + 

K D1 P 

bBb 

+ K B2 P 2 

+ K b1 P 3 

* K B4 P 4 

+ K B5 P 5 

+ K^P 
+ K B4 P 4 

* K B5 P S 

♦ K B6 P 6 

* K B7 3p 7 

+ K D4 p 
+ d6 P 6 

+ <B7 P 7 

. K D4 p 
+ ^B8 P 8 

* K^Pg 

+ K P ^P 
+ ^B8 8 

+ k B9 P 9 

+ K P ^P 
+ ^BB P B 



'D2 p 
W B 

'D3 p 
W B 

,D4 p 

'BB r B 


( 11 ) 


co C D1 f D2 .03 -04 r RS 

S2 = S B + S B + S B + 6 B + S B 


( 12 ) 


S3 = S 


TB 

B 


(13) 


AB - k Jb p b 


(14) 


Rearranging equation (11) gives 


ci l^DI d [.Dl .D2 i D I.Dl |^D2 \ D 
bi = + f' B2 + ^b 2 / p 2 + PB3 + ^B3 | P 3 


.D2 J3 D .D2 ,D3 p .D3 .D4 L 

+ K B4 + K B4 P 4 + K B5 + K B5 P 5 + K B6 + K B6 ) P 6 


.D3 J4 L .D4 ,RS L .04 .RS D 
+ B7 + ^B7 / 7 + B8 + B8 P 8 + «9 + ^B9 P 9 


.Dl .D2 .D3 .04 .RS D 

+ BB + ^BB + ^BB + BB + ^BB B 


(15) 
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and equation (10) is rewritten as 


SI + AB = AL - S2 + S3 


(16) 


and substitution of equations (12) to (15) into equation (16) yields the 10th 
compatibility equation which is 

„D1 d [ v D1 „D2Y („ Dl D2\ („D2 

K B1 P 1 + K B2 + K B2 ] P 2 K B3 + K B3 P 3 + K B4 + K B4 P 4 


yD2 „D3 D LD3 ^04 L L,D3 .D4 L 

+ K B5 + K B5 P 5 + K B6 + K B6 P 6 + K B7 + K B7 P 7 


f 


D4 


,RS 


,04 


,RS 


,D1 


,D2 


K B8 + K B8 P 8 + K B9 + K B9 P 9 + K BB + K BB + K BB T "BB T ^BB /' B 


,D3 „D4 ,,RS 1 D 

^nn + ^nn Jr r 


J B D - A I * D1 ^ D2 * 03 <D4 x RS JB 

+ K BB P B " AL S B + S B + S B + S B + S B + S B 


(17) 


Each of the components that make up the compressor rotor assembly were 
separately analyzed with the FEATS program in order to determine the j and 
the Si terms in the compatibility equations. The components were first ana- 
lyzed under centrifugal and thermal loading conditions only (effects of the 
forces. Pi, were neglected) in order to determine the Si terms. The blade 
loads were applied as distributed nodal point loads and the temperature distri- 
bution over the components was manually input by prescribing a temperature for 
each element. These elemental temperatures were obtained from a previous 
steady-state heat transfer analysis. The Si terms for the compatibility 
equations were then determined by noting the resulting nodal displacements at 
the contacting surface(s) of each component. This procedure required a total 
of seven separate computer runs. Next, the k^j terms (influence coeffi- 
cients) were determined by loading each component with a trial load at a con- 
tact surface (j) and noting the resulting nodal displacement at some contact 
surface (i). No other loads were applied to the component. Several trial 
load conditions had to be applied to each component in order to obtain all of 
the influence coefficients. This was the most time-consuming portion of the 
analysis since it required 24 separate computer runs to determine the 102 
influence coefficients. 


Since the pilot interfaces in the assembly are designed with a 0.0254-mm 
(0.001-in.) radial interference fit, the Ay^ terms in the compatibility equa- 
tions were set to 0.0254 mm (0.001 in.). Notice that this is a differential 
deflection term which can be placed on either side of the equation according 
to a chosen sign convention. 

The desired tie bolt preload was chosen to be 47 148 N (10 600 lb). For 
the given bolt size, this load results in a bolt elongation of approximately 
0.508 mm (0.020 in.). Therefore the AL term in the 10th compatibility equa- 
tion (eq. (10)) was set to 0.508 mm (0.020 in.). 

Once the general formation of the 10 compatibility equations was accom- 
plished, the equations were rearranged in matrix form by transforming unknown 
terms, P ^ , and their coefficients to the left side, and the known terms to the 
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right side. An initial solution of the set of equations yielded negative val- 
ues for some of the axial forces; therefore those forces were set equal to zero 
and another solution was obtained. After several iterations, it was found that 
the forces ? 2 , P 4 , P 6 » and Pg were all equal to zero. The values of the 
forces in the radial direction all turned out to be positive, which indicates 
that no gaps will develop at the pilot interfaces during operation. Also, the 
load Pg was positive, which means that the tie bolt remains loaded during 
operation. 

After the interacting forces on the assembly were determined, another 
individual finite element analysis of the components was performed, including 
these forces, to determine the final operating stresses of the compressor 
rotor. 


IMPROVED ANALYSIS METHOD 

For this analysis technique, a two-dimensional finite element model of 
the entire compressor was developed which incorporates all of the components 
and also properly simulates the effects of all contacting surfaces. The analy- 
sis, which uses the MARC finite element computer program (ref. 8 ), was done in 
two parts. The first part was a heat transfer analysis that determined the 
steady-state temperature distribution over the compressor due to the airstream 
temperatures and bearing operating temperatures. The resulting elemental tem- 
peratures were saved on a computer tape and later incorporated into the com- 
puter analysis that was performed for the stress-compatibility analysis. 

The finite element mesh shown in figure 5 is comprised of 371 eight-node 
quadrilateral axisymmetric elements. The mesh actually consists of six sepa- 
rate finite element models joined together by MARC gap elements; the gap 
elements were placed at the surfaces that are in contact. There are 10 con- 
tacting surfaces in the model, 5 are the pilot interfaces, and 5 are the sur- 
faces that are held tightly together by the force of the tie bolt. These 10 
surfaces were modeled with 20 fixed-direction gap elements, 2 gap elements per 
surface. Ten of the gap elements were orientated in the radial direction and 
model the pilot contact surfaces, while the other 10 were orientated in the 
axial direction and model the surfaces that are held together by the tie bolt. 
Note that the tie bolt was integrated within the rear shaft since the assump- 
tion was made that the threaded connection of the tie bolt is rigid. 

The gap elements can be defined with either a closure distance or an 
interference fit. If a closure distance is chosen, the contacting surfaces 
are initially assumed to be separated by that distance (i.e., a gap exists). 

On the other hand, if an interference fit is chosen, the surfaces are initially 
assumed to be in contact. In this case, the MARC program simulates the pres- 
sure loading on the surfaces by forcing the two end nodes of the gap element 
(there is one node on each surface) to be displaced from each other an amount 
equal to the interference fit. As a result, interference fits on the model 
cause a stress distribution before any of the operating loading conditions are 
applied. In this analysis, only interference fits were used for the gap 
elements . 

The radial interference fit at each of the pilots was set by specifying a 
0.0254-mm (0.001-in.) interference in the radial direction gap elements. The 
preload of the tie bolt was set by specifying an interference fit in all of the 
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axial gap elements. The desired tie bolt preload was chosen to be 47 148 N 
(10 600 lb); for the given bolt size, this load results in a bolt elongation 
of approximately 0.508 mm (0.020 in.). Since there are five axial contact sur- 
faces, each of the gap elements at these surfaces were given an interference 
fit of 0.1016 mm (0.004 in.). This method seems to represent the effects of 
the preloaded tie bolt quite accurately. 

As boundary conditions, all of the nodes along the axis of rotation were 
constrained against radial motion. Also, the node points located at the con- 
necting bolt diameter, between the front shaft and the first disk, were tied 
together in the axial direction. This simulates a rigid bolted connection. 

The stress-compatibility analysis of this system requires a nonlinear 
solution, therefore the system equations must be solved incrementally by the 
MARC finite element program. The governing equations are expressed in an 
incremental form as 


K du = df 

where du and df are incremental displacement and force vectors respectively 
and K is the stiffness matrix. The model was loaded with nine successive 
equal load increments. After the ninth load increment was applied, the total 
load on the model was equivalent to the actual loads at 47 500 rpm. The load- 
ing increments were specified within the input file such that the finite ele- 
ment program applies them automatically in one computer run. The element 
temperatures were applied to the model before the first load increment by read- 
ing them off of the computer tape created in the heat transfer analysis. One 
ninth of the centrifugal load and blade loads was then applied at each of the 
following nine load increments. The blade loads were applied as distributed 
nodal point loads. 

If a gap develops during one of the load increments, the stiffness matrix 
changes and must be recalculated. The MARC program uses the full Newton- 
Raphson method as the solution procedure to determine the new stiffness matrix. 
Within each load increment, the program performs iterations, recalculating and 
refactorizing the stiffness matrix until the convergence criterion is met. 
Convergence is met when the ratio of the maximum residual force to maximum 
reaction force is less than a given tolerance. 

The complete analysis required only two computer runs: one for the heat 

transfer analysis and the other for the stress-compatibility analysis. The 
program calculated the nodal displacements and gap element conditions after 
every load increment. The element stresses were calculated after the ninth 
(final) increment only, although they could have been asked for at any or all 
of the increments. 


RESULTS AND DISCUSSION 

The flexibility analysis method utilized the FEATS finite element compu- 
ter program in order to obtain certain values within the set of compatibility 
equations and also to perform the final stress analysis of the rotor compo- 
nents. This effort required a total of 31 separate finite element computer 
runs. A simultaneous solution of the matrix of compatibility equations had to 
be made three times to solve for the final unknown forces. 
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The improved analysis method, or gap element method, required only one 
finite element computer solution to determine the operating stresses and to 
determine the compatibility of all contacting surfaces for the compressor 
rotor. There were no compatibility equations to solve. A steady-state heat 
transfer analysis was done prior to the stress analysis in order to determine 
the element temperatures that were used as the thermal loading conditions. 

This heat transfer analysis did not require much time or effort since the same 
finite element mesh was used as in the stress-compatibility analysis. 

In both analysis methods, during operating conditions, the forces P£, 

P 4 , P 5 , and Pg were all found to be equal to zero, which indicates that there 
are no moment forces between components. The pilot surfaces were all deter- 
mined to be closed since the forces P-j , P 3 , P 5 , P 7 , and Pg all had positive 
values. Also, the tie bolt was tight since Pg had a positive value. The 
forces resulting from each analysis method are shown in table II. The values 
given for these forces are integrated around the circumference of the pilot 
contact surface. The tie bolt load, Pg, decreased from 47 148 N (10 600 lb) 
to approximately 33 800 N (7600 lb) during operation. This was expected since 
the Poisson's ratio effect slightly reduced the axial dimension of the assem- 
bly of components. The forces determined by the flexibility method compare 
quite well to those determined by the improved method with the largest discrep- 
ancy seen in the value for P]. One explanation for the differences may be 
that the improved method used an eight-noded element, which is much more accu- 
rate than the four-noded element that was used in the other method. Also, 
there are differences in the solution algorithms of the two finite element 
codes that were used. 

Contour plots of the radial and tangential stress distribution on the com- 
ponents are shown in figures 6 to 11 for the flexibility method and in figures 
12 to 17 for the improved method. The stress contours resulting from each 
analysis method compare very well to each other. The effects of the radial 
pilot loads can be seen by noticing that the radial stress is compressive (neg- 
ative) at the pilot contact surfaces. Also, from examination of the elemental 
stress values, it was found that the radial stress is continuous across the 
pilot interfaces at the point of contact, which is a required boundary condi- 
tion (ref. 9). Generally, the peak tangential stress values on the disks are 
higher using the improved method. This can be attributed to the use of the 
eight-noded elements, which give better results in areas that have high stress 
gradients. 


CONCLUDING REMARKS 

The stress-compatibility results that were obtained by using the improved 
analysis method agree very closely with those that were obtained by using the 
flexibility analysis method. In both cases, the design of the compressor 
rotor is satisfactory for operation within the stated operating conditions. 

The disadvantage of the flexibility analysis method is that it is more 
time consuming and tedious than the improved analysis method. Many finite ele- 
ment solutions had to be performed to determine data that was needed in the 
compatibility equations. This analysis method, being so tedious, is usually 
only performed at the maximum operating speed. Therefore, information such as 
displacements or pilot movement trends cannot be determined at lower speeds 
unless another analysis is performed. If the analysis determines that any or 
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all of the pilot contacts have opened during compressor operation, then the 
pilots have to be redesigned, possibly by reversing their orientation. This 
requires the formulation of new finite element models representing the modi- 
fied components, and the complete analysis has to be performed again. In some 
cases, this cycle had to be completed four or five times before a satisfactory 
design was di scovered. 

In the case of the improved analysis method, redesigns are accomplished 
with relative ease, requiring only minor modification of the finite element 
model and a solution by the MARC program. Another advantage is that the force 
at each gap element can be output at the end of each load increment. There- 
fore, as the compressor rotor speed increases, the engineer can determine if 
the pilot contact surfaces are tending to tighten, which is a desirable condi- 
tion. The flexibility analysis method only reveals the contact forces at the 
speed that the analysis was performed. 

Given the ease, timeliness, and reliability of the improved analysis 
method, it is a preferred alternative when analyzing multicomponent rotating 
systems . 
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TABLE I. - MATERIAL PROPERTIES OE COMPRESSOR COMPONENTS 3 


A1 1 oy 

Temperature 

Modul 

is of 
ci ty 

Poisson's 

ratio 

Coefficient of 
thermal expansion 

°C 

°F 


N/m 

psi 

m/m °C 

in. /in. °F 

T i -6A1 -4V 

-18 

204 

0 

400 

5.75x10^ 

5.40x10’ 

16.5xl0f 

15.5x10° 

0.33 

.33 

8.82xl0-f 

9.36x10-° 

4.90x10-5 

5.20xl0- 6 

AISI 4340 
alloy steel 

-18 

204 

0 

400 

10.45x10^ 

9.41x10’ 

30 .OxlOf 
27.0x10° 

0.30 

.30 

11.16x10-5 

13.32x10-5 

6.20xl0- 6 

7.40x10-6 

18-200 maraging 
steel 

-18 

204 

0 

400 

9 . 13x1 0^ 
9. 13x10’ 

26.2x10^ 

26.2x10° 

0.264 

.264 

10.1x10-5 

lO.lxlO -6 

5.60x10-5 

5.60xl0- 6 

1 

A1 1 oy 

Thermal conductivity 

Mass density 


W/m °C 

Btu/sec-in. °F 

kg/m^ 

lb sec^/in.^ 

T i -6A1 -4V 

0.806 

1.027 

9.26x10-5 
1 . 18x1 0 -4 

4423 

4423 

4.14xl0“4 

4.14x10-4 

AISI 4340 
alloy steel 

4.03 

4.63 

4.63x10-4 

5.32xl0- 4 

7831 

7831 

7.33x10-4 

7.33x10-4 

18-200 maraging 
steel 

4.35 

4.35 

5.00x10-4 

5.00x10-4 

8002 

8002 

7.49x10-4 

7.49x10-4 


a From reference 10. Properties are linearly interpolated between the temperature limits. 


TABLE II. - PILOT SURFACE FORCES PREDICTED BY TWO 
ANALYSIS METHODS 


Load 

Force predicted by 
fl exi bi 1 i ty 
analysis method 

Force predicted by 
improved 
analysis method 


N 

lb 

N 

lb 

fj 

43 928 

9 876 

29 392 

6 608 

67 756 

15 233 

64 215 

14 437 

p 5 

77 360 

17 392 

76 857 

17 279 

P 7 

67 614 

15 201 

67 432 

15 160 

p 9 

149 221 

33 548 

139 520 

31 367 

P B 

34 107 

7 668 

33 813 

7 602 


AXIAL CONTACT 



FIGURE 1. - TYPICAL PILOT DESIGN USED IN ROTATING MACHINERY. 
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DISKS ^ 








6 


CONTOUR 


RADIAL 

STRESS, 



(A) RADIAL STRESS. 




MPa 

PSI 

A 

350.5 

50 saa 

B 

321.6 

A6 655 

C 

292.7 

A 2 A65 

D 

263.9 

38 275 

E 

235.0 

3A 086 

F 

206.1 

29 896 

G 

177.2 

25 706 

H 

1A8.3 

21 517 

I 

119. A 

17 327 

J 

90.6 

13 137 

K 

61.7 

8 9A8 

L 

32.8 

A 758 

M 

3.9 

569 

N 

-2A.9 

-3 620 


CONTOUR 

TANGENTIAL 


STRESS, 


MPa 

PSI 

A 

361.2 

52 395 

B 

3A7.3 

50 373 

C 

333.3 

A8 351 

D 

319. A 

A6 329 

E 

305. A 

AA 307 

F 

291.5 

A2 285 

G 

277.6 

AO 263 

H 

263.6 

38 2A1 

I 

2A9.7 

36 219 

J 

235.7 

3A 197 

K 

221.8 

32 175 

L 

207.9 

30 153 

M 

193.9 

28 131 

N 

180.0 

26 109 

0 

166.0 

2A 087 

P 

152.1 

22 065 

Q 

138.2 

20 0A3 


(B) TANGENTIAL STRESS. 

FIGURE 8. - STRESS CONTOURS FOR SECOND DISK PREDICTED BY FLEXIBILITY ANALYSIS METHOD. 




'—A 


(B) TANGENTIAL STRESS. 


CONTOUR RADIAL 

STRESS. 



MPa 

PS I 

A 

389.2 

56 959 

B 

361.9 

52 926 

C 

333.6 

98 397 

D 

305.9 

99 369 

E 

278.1 

90 391 

F 

250.3 

36 312 

G 

222.6 

32 289 

H 

199.8 

28 255 

I 

167.0 

29 227 

J 

139.2 

20 198 

K 

111.5 

16 170 

L 

83.7 

12 192 

M 

55.9 

8 113 

N 

28.2 

9 085 

0 

.9 

56 

P 

-27.9 

-3 971 


CONTOUR TANGENTIAL 

STRESS. 



MPa 

PSI 

A 

389.0 

56 925 

B 

375.0 

59 900 

C 

361.1 

52 375 

D 

397.1 

50 350 

E 

333.1 

98 325 

F 

319.2 

96 300 

G 

305.2 

99 275 

H 

291.3 

92 250 

I 

277.3 

90 225 

J 

263.3 

38 200 

K 

299.9 

36 175 

L 

235.9 

39 150 

M 

221.5 

32 125 

N 

207.5 

30 100 

0 

193.5 

28 075 

P 

179,6 

26 050 


FIGURE 9. - STRESS CONTOURS FOR THIRD DISK PREDICTED BY FLEXIBILITY ANALYSIS METHOD. 
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CONTOUR RADIAL 

STRESS, 



MPa 

PSI 

A 

<116.2 

60 366 

B 

388.4 

56 345 

C 

360.7 

52 325 

D 

333.0 

48 305 

E 

305.3 

44 284 

F 

277.6 

40 264 

G 

249.9 

36 244 

H 

222.1 

32 223 

1 

194.4 

28 203 

J 

166.7 

24 183 

K 

139.0 

20 162 

L 

111.3 

16 142 

M 

83.6 

12 122 

N 

55.8 

8 101 

0 

28.1 

4 081 

P 

.4 

61 

Q 

-27.3 

-3 959 


CONTOUR TANGENTIAL 

STRESS, 



MPa 

PSI 

A 

383.1 

55 571 

B 

365.6 

53 037 

C 

348.2 

50 504 

D 

330.7 

47 970 

E 

313.2 

45 437 

F 

295.8 

42 903 

G 

278.3 

40 369 

H 

260.8 

37 836 

I 

243.4 

35 302 

J 

225.9 

32 768 

K 

208.4 

30 235 

L 

191.0 

27 701 

M 

173.5 

25 168 

N 

156.0 

22 634 

0 

138.6 

20 100 

P 

121.1 

17 567 


(B) TANGENTIAL STRESS. 

FIGURE 10. - STRESS CONTOURS FOR FOURTH DISK PREDICED BY FLEXIBILITY ANALYSIS METHOD. 
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